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1 - Introduction 

The theory of weak solutions for nonlinear conservation laws is now well developed in the case 
of scalar equations [3] and for one-dimensional hyperbolic systems [1, 2]. For systems in several 
space dimensions, however, even the global existence of solutions to the Cauchy problem remains 
a challenging open question. In this note we construct a conterexample showing that, even for a 
simple class of hyperbolic systems, in two space dimensions the Cauchy problem can be ill posed. 

The systems on that we consider take the special form 

+ =° i = l,...,n. (1.1) 

Q = l " 

For physical motivations, see [2] or [5]. Given a sufficiently regular initial condition 

u[Q,x) = u{x) , (1.2) 
the solution of the Cauchy problem can be constructed as follows. 

1. Set p = \u\ and solve the Cauchy problem for a scalar conservation law 

+ =0' p{Q,x) = \u{x)\. (1.3) 

This will provide the absolute value of the solution. 

2. To find the angular component 9 = u/|u|, consider the O.D.E. 

± = f{p{t,x)) (1.4) 

with / = (/i, . . . ,fm)- Assuming that (1.4) has a unique solution for almost every initial data 
x(0) = y, denote hy 1 1-^ x{t) = ^tU the corresponding trajectory. 

3. Given t > and x G iR*", call y = <I>~*x be the (unique) point such that ^tU = x. Assuming 
that the measurable maps are well defined (up to sets of zero measure), the angular component 
of the solution of (1.1)-(1.2) is then obtained as 

e{t,x) = e{^-^x) , (1.5) 



where 6 = u/\u 



If the solution p = p{t,x) of the scalar conservation law (1.3) remains smooth (or piecewise 
smooth), it is easy to check that the above procedure actually yields a weak solution u to the Cauchy 
problem (1.1)-(1.2). In general however, even if the initial data p{0,-) is smooth, regularity may 
be lost in finite time. It is thus interesting to understand what kind of assumptions can provide 
the existence of entropy weak solutions. A natural set of conditions is the following: 

(Al) The flux function F(p) = f{p) p is Lipschitz continuous. 

(A2) The initial data is bounded, measurable, and bounded away from zero, i.e. 

uGL°°, < a < [^(a;)! < 6 for a.e. a; G iR"* . 

In the case of one space dimension, the above assumptions yield the global existence of a 
unique entropy weak solution to the Cauchy problem. Indeed, by the fundamental theorem of 
Kruzhkov [3], the equation (1.3) admits a unique entropy weak solution p = p{t,x), satisfying 

a < p{t,x) < b. (1.6) 

Let Xi{t) < X2{t) be any two solutions of (1.4). For every t > the conservation equation implies 




By (1.6) this yields 

^(x2(0)-a;i(0)) < X2{t)-xi{t) < ^(^2(0) - a;i(0)). 

Hence the maps 

{t,y)<-^^ty, {t,x)<-^^-*x 
are both well defined and Lipschitz continuous. It is now clear that (1.5) yields the desired solution. 

The example constructed below will show that, under exactly the same assumptions (Al)- 
(A2) which provide the existence and uniqueness of solutions in the one-dimensional case, in two 
space dimensions the Cauchy problem (1.1)-(1.2) is not well posed. The reason is that, although the 
scalar conservation law (1.3) always admits a unique entropy solution, the corresponding O.D.E. in 
(1.4) has discontinuous right hand side. For a given time t > 0, the flux function y ^ty may 
not admit a well deflned inverse, because of oscillation phenomena. 



2 - A counterexample 

(0,1) the canonical basis on M^. 
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Observe that 
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We can write F{p) = f{p) p, with 
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In particular, 



/(2) = -2^1 



/(3) = 0. 



/(4) = . 



We now look at the O.D.E. 



X = f{p{t,x)) 

in connection with special types of solutions p = p{t, x) of the conservation law 

Pt + V, •F(p) = 0. 

For a given set Q C M^, we consider three entropy weak solutions of (2.6), namely 

p\t,x) = 3, 



p'{t,x) = 
p\t,x) = { 



4 ii x-te2 eQ, 

3 otherwise, 

2 if x - tei G Q , 

3 otherwise. 



figure 1 



We begin by studying what happens to the corresponding trajectories of (2.5) in the special 
case where Q is a rectangle: 

Q = [0, a] X [0, b]clR^. 

1. In the case p = p'', we trivially have 

x = ^\y = y. 

2. When p = p^, we have (fig. 1) 

^^|e2/4 if a;iG[0,a], X2-te[0,b], 

1 otherwise. 

Points on the vertical strip 

Ta = {(xi,X2) ; xi G [0,a]} 

are eventually shifted to the right by an amount bei/3, while all other points do not move. More 
precisely: 

$«y-y = Pei/3 if yie[0,a], y2>b, t>iy2-b)+4b/3, 

1 otherwise. 

This is illustrated in fig. 1. One has p = 4 on Q and p = 3 outside. The rectangle Q moves 
upward with unit speed. At any given time t, points which lie inside the rectangle move with 
speed X = 62/4. The points which initially lie on the segment 7 are eventually displaced onto the 
set $"7. 

3. When p = p\ we have 

i; = /-ei/2 if xi-tG[0,a], X2G[0,6], 

1 otherwise. 

Points on the horizontal strip 

Pfe = {{xi,X2) ; X2 G [0,6]} 

are eventually shifted to the left by an amount — aei/3, while all other points do not move. More 
precisely: 

^\y -y=\ if y2 e [0> b], yi>a, t> {yi - a) + 2a/3 , 

1 otherwise. 

Starting with these elementary solutions, we can construct more complex ones. For example, 
in (2.7) can take 

Q = Q" = U [k,k + l] X [k,k + 3] 

k even 

In the corresponding equation (2.5), this generates the shift 



1^2 if [[yi]] is even, 

^ ^ ^ ~ \ if [[yi]] is odd. 



^1 

figure 2 



On the other hand, in (2.8) we can take 

k even 

In the corresponding equation (2.5), this generates the shift 

^K, _,,-[ if [[2/2]] is even, 

^ ^"10 if Wis odd, 

where [[s]] denotes the largest integer < s. 

In the following, we write the coordinates of a point x = {xi,X2) € in binary digits, say 

Xl = 5^afe2-^ X2 = 5;/3fe2-^ (2.9) 

k k 

with ak, Pk G {0, 1}. For every A;, the previous analysis shows that one can construct sets Q\., 
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QI. contained inside disjoint strips 

QIq{{xi,X2); X2+X1 g2-'=[28, 32]}, 



Qi^{{xi,X2}; X2 + X1 G2-'=[22, 27]}, (2.10) 



Ql^{{xi,X2); + G2-'=[16, 21]}, 
and such that the following holds, 
(i) Given the solution 

pUt,x) = \t if ^-*e2GQl, 

''^ ^ [3 otherwise, 

the corresponding equation (2.5) generates the shift 

^K-v = ( ^"'^2 if «fc = 0, 

^ ^ lo if ak = l. 



(ii) Given the solution 



Pkit^x) = 

the corresponding equation (2.5) generates the shift 



2 if X + tei G , 



k 

3 otherwise. 



-2-^ei if Pk + 

if = 



(iii) Given the solution 

pi(i,:r) = |4 if x-te2GQ«, 

^ i 3 otherwise, 

the corresponding equation (2.5) generates the shift 



- y = 2-'=-^e2 if = 0, 

if afc = l. 



As shown in fig. 2, the set QI. consists of an array of (gray) rectangles of size 2 x 3 ■ 2 
the set Q\. consists of an array of (black) rectangles of size 2^"^^ x 2~^, while the set consists of 
an array of (gray) rectangles of size 2"*^ x 3 • 2^"*^. Moving the gray rectangles upward with unit 
speed and moving the black ones to the right with unit speed, the three arrays will never overlap. 
We can thus consider the composite solution 

'4 if X - te2 G Qi U Q« , 
Pk{t,x)=i2 if x-tei^Ql, (2.11) 

, 3 otherwise. 
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figure 3 



The differential equation 

X = f{pk{t,x)) 

will eventually shift points on the plane, according to the composition 

*, = $«o$^o$«. (2.12) 

This is illustrated in fig. 3. Each square here has sides of length 2~*^~^. In all four pictures, 
the asterisque marks the samepoint in the plane. Performing first the upward shift , then the 

leftward shift and finally the upward shift <1>^ , one obtains a measure-preserving transformation 
whose main property is the following. In connection with the dyadic decomposition (2.9), define 

Xj = {{xi,X2); 0j = O}, X'j = {{xi,X2); 0j = l}. 

Then, up to sets of measure zero, one has 



^k{Xk) = Xk+i, *fc(X^) = X^+i. (2.13) 

Intuitively, one can think at the moving sets Q^, Q\., Q\ as a comb, moving in the north-east 
direction. The effect of its passage is to displace points on the set Xk onto the set Xk+i- Putting 
an arbitrarily large number of these combs one next to the other (fig. 4), we now construct a 
sequence of solutions where the initial data converge strongly in Lj*^^ but at a later time t the 
corresponding solutions have no strong limit. Hence the limiting Cauchy problem cannot be well 
posed. 
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Define a solution of (1.3) by setting 

f4 if x-ie2GUo<fc<n(Q[ugi) 

2 if a;-tei GUo<fe<nQfc' 

^ 3 otherwise. 



Pn{t,x) = < 



(2.14) 



Consider the Cauchy problem (1.1)-(1.2) with / given at (2.3). We choose an initial value CZ„ of 
the form 

Un{x) = p„(0,x) • (cos^(a;) , sin^(a;)) (2-15) 



with 



e{x) = 



(3 if [[a;2]] is even, 

—(3 if [[2:2]] is odd, 



for a fixed angle /3 G ]0, 7r/2[ . We claim that the sequence of weak solutions has no limit in Lj^^ 
as n — > 00. Indeed, for 32 < xi + X2 < ^ one has 



Un{t,x) = 3{cos9n{t,x) , siu 6'„(t, x)) , 



with X = (xi,X2) and 



i{t,x) 



13 if [[2'"X2]] is even, 

-/3 if [[2"a;2]] is odd. 



As a result, when n ^ 00, the Un do not converge in Lj^^^ and the Cauchy problem is not well 
posed. 




figure 4 



On the other hand, still in the range where 32 < xi + X2 < t, one can compute the weak limit 
of the sequence {un}, given by 

u{t,x) = (3COS/3, 0). 
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To check whether this weak Umit provides a distributions solution to the system (1.1), (2.3), it 
suffices to examine what happens in a neighborhood of the hne xi + X2 = t. Computing the time 
derivative 




for a set of the form 

Q{t) = {{xi,X2) , X2e[a,b], xi G [t-X2-2-^,t-X2 + 2-'=]}, 

one checks that the conservation equations are satisfied provided that F(3cos/3) = F{3) = 0. In 
other words: 

- If cos (3 < 1/3, then the weak Hmit u provides a weak solution to the system of conservation laws. 
This yields a curious example where entropy is dissipated by a linearly degenerate field. 

- If cos/3 > 1/3, then the weak limit u is not weak solution. In this case, it is not even clear 
whether any solution exists at all. 



3 - Concluding remarks 

1. All systems of the form (1.1) are hyperbolic, but in a certain way pathological. In one space 
dimension, the system 

nt + (/(M)u)^ = (3.1) 

takes the quasilinear form 



ut + A{u)u^ = 0, A{u) = f{\u\) I + f'{\u\) 



At a given point n 7^ 0, the Jacobian matrix A{u) has the 1-dimensional eigenspacc E = Mu 
spanned by the vector u, corresponding to the eigenvalue A = /(|n|) + \u\. Moreover, 

one finds a second eigenvalue A* = of multiplicity n — 1, corresponding to the orthogonal 

eigenspace E-^ (fig. 5). 

Even in one space dimension the Riemann problem can have multiple solutions, all obtained 
as vanishing viscosity limits. For example, consider 

Ut + {\u'^\u)^ = ueM'^, (3.2) 

^/Q x) = h^ if X < 0, 

An entropic solution is provided by 

u{t, x) = p{t, x)ei , 
where p is the entropy weak solution of the corresponding scalar problem 
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figure 5 



This is clearly a limit of vanishing viscosity approximations [3] . 
In addition, there is the second admissible solution 

u{t,x) = h^ ^ (3.3) 

^ ' ^ \ -ei if x>t. ^ ' 

Indeed, consider any smooth function 6 : M i-^ [0,Tr] such that 

[ vr if X > 1. 

Then the Cauchy problem for (3.2) with initial data 

Un{0,x) = (^cos6{nx) , sm9{nx)^ 

has a unique smooth solution, namely Un{t,x) = u{0,x — t). This is also a limit of vanishing 
viscosity approximations. Letting n — > cx), these solutions converge to (3.3). 

We remark, however, that our counterexample is not related to the presence of a singularity 
at the origin. Indeed, the initial data at (2.14)-(2.15) are contained in a bounded domain whose 
convex closure does not contain the origin. 



2. The flux function is not smooth but only Lipschitz continuous. The oscillations in the solution 
of the scalar conservation law (1.3) are not damped, because of the linear degeneracy of the flux. 
This indeed leaves open the possibility that, with a smooth flux, such a pathological behavior will 

not be observed. 

We recall that, in the one dimensional case, it is common to replace a general flux by a 
polygonal approximation F^,. This is indeed the basis for the method of front-tracking approxima- 
tions. In one space dimension, the Cauchy problems related to all polygonal fluxes are well posed 
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[1]. Taking a limit, one thus obtains a proof of well posedness for the original Cauchy problem. 
The multidimensional case is thus completely different. 



3. What our example points out is that, to get well posedness, the initial data should be chosen 
in a space of functions smaller than L°°. In view of the one-dimensional theory, it is natural to 
conjecture that a global weak solution should exist for initial data u with bounded variation. A 
proof of this conjecture (far from easy!) could rely on a compactness property for fluxes generated 
by O.D.E's with right hand side in BV. More precisely, consider a sequence of smooth maps 
: [0, T]xlR^ ^ IRT such that 

\U{t,x)\<Cr, (3.4) 



BV 



T 
J 



i=l 



_d_ 

dx, 



dxdt < C2 ■ 



(3.5) 



Call t x{t) = (y) the solution of 

X fvi^i 23), 



x{Q) = y. 



Moreover, assume that the fluxes are all nearly incompressible, so that, for every bounded set 
A c iR™, 

-^meas(A) < meas ($4'' (A)) < C3meas(yl), (3.6) 
for some constant C3 and alH G [0, T], 1/ > 1. 

Conjecture: By possibly exstracting a subsequence, one has the convergence 

in Li„, 

for some measurable flux also satisfying (3.6). 

Notice that, by a well known compactness theorem, for some / G L°° we certainly have 



in Moc 



The existence of a limit flux $ would provide a new existence result concerning the Cauchy problem 
for the discontinuous O.D.E. 

x = f{t,x), x(0)=y (3.7) 



valid for a.e. initial data y G IRJ^, somewhat extending the result in [3]. 
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